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We study the fermionic spectral function in a holographic superconductor model. At zero tem- 
perature, the black hole has zero horizon and hence the entropy of the system is zero after the 
back reaction of the condensate is taken into account. We find the system exhibits the famous 
peak-dip-hump lineshape with a sharp low-energy peak followed by a dip then a hump at higher 
energies. This feature is widely observed in the spectrum of several high-T c superconductors. We 
also find a linear relation between the gap in the fermionic spectrum and the condensate, indicating 
the condensate is formed by fermion pairing. 



Introduction-The proper description of the physics of 
strongly interacting fermions has long been the major 
issue for our understanding of the quark-gluon interac- 
tions and condensed matter systems near quantum crit- 
icality. Recent development of the holographic corre- 
spondence between a gravitational theory and a quan- 
tum field theory, which first emerged under the anti- 
de-Sitter space/conformal field theory (AdS/CFT) cor- 
respondence P, [2, El , has provided a useful framework 
to describe systems with strong interactions (see e.g. 

BHiMMIIII). 

Recently, a gravitational model of hairy black holes 
[H, EH ! where the Abelian symmetry of Higgs is spon- 
taneously broken below some critical temperature, has 
been used to model high-T c superconductivity (HTSC) 
[Til ]. On the fermion side, it has been shown that 
Fermi liquids emerge by deforming away from the quan- 
tum critical point by increasing the fermion density [l5| . 
Also, a theory describing non-Fermi liquids, which de- 
scribes the normal state of the high-T c superconducting 
cuperates and metals near quantum critical points, has 
been proposed 16]. The response functions of compos- 
ite fermionic operators in a class of strongly interacting 
quantum field theories at finite de nsity are computed us- 
ing the AdS/CFT correspondence [H,[I3- In particular, 
in Ref. [l6[ the Dirac equation is probed in the back- 
ground of an extremal AdS4 charged black hole which 
possesses a finite horizon area and thus finite entropy 
at temperature T = 0. Analysis of the fermionic spec- 
tral function gives gapless fermionic excitations at dis- 
crete shells in momentum space, signaling the existence 
of a Fermi surface, and the spectral weight exhibits novel 
properties, for example, particle-hole asymmetry. 

In this letter, we study the fermionic excitations of a 
holographic superconductor. This is motivated by the 
great importance to understand how fermionic quasipar- 
ticle excitation is affected by the condensate, as this may 
hold the key to the mystery of HTSC. As an intrinsic 
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strong coupling system, HTSC has posted great theoret- 
ical difficulties since its discovery. With the AdS/CFT 
correspondence, where a strongly coupled field theory can 
be mapped to a weakly curved gravity theory, nonpertur- 
bative computations in quantum field theory turns into 
solving classical Einstein gravity. At T = 0, the sym- 
metry broken phase has degenerate vaccua. However, 
since only one state is accessible, it is more reasonable to 
choose a gravitational model with zero horizon at T = 0. 
This requirement can be achieved by including the full 
back reaction [l8[ of the condensate to the gravity geom- 
etry [l9j . By applying the fermion spectral function tech- 
nique developed in [201 ] , we find that the system exhibits 
the famous peak-dip-hump (PDH) lineshape with a sharp 
low-energy peak followed by a dip then a hump at higher 
energies. This feature is widely observed in the spec- 
trum of several HTSC, for example, in the angle resolved 
photoemission spectroscopy (ARPES) data near (n, 0) 
point of the Brilluoin zone in Bi2Sr2CaCu 2 Os(Bi2212) 
[2lL [22)], and in tunneling experiments [231 ] . This fea- 
ture has been attributed to the many-body effects of the 
fermionic quasiparticle coupling to a bosonic collective 
mode, such as the (ir, ir) resonance mode in inelastic neu- 
tron scattering [24], and it is argued to be related to 
the microscopic pairing interaction. Understanding the 
physics of this PDH lineshape feature is therefore very 
important for the construction of a microscopic mecha- 
nism for HTSC. Furthermore, by assigning the fermion 
charge half that of the condensate, we find a linear rela- 
tion between the gap in the fermionic spectrum and the 
condensate. This result indicates that the condensate is 
related to the fermion pair in our theory. Since we use 
fermionic degrees of freedom as a probe to the conden- 
sate, this might provide us a pairing mechanism without 
a "glue" for superconductivity. 

Gravitational Model of the Superconductor-Inspired by 
the AdS/CFT correspondence, we conjecture that there 
is a strongly interacting CFT dual to the gravitational 
model that we will describe below. It is assumed that 
even though we do not know what exactly the underlying 
field theory is, we can still study its properties from the 
gravitational side, using the AdS/CFT machinery. 

We start with the following action with a vector field 
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where the gravitational coupling k 2 = 8itGn, m and q 
are the mass and charge of the scalar field. We will use 
the units in which the radius of curvature L = 1. A fully 
back-reacted ansatz for the metric is 



dr 2 



ds 2 = -g(r)e- x(r) dt 2 + ^ + r 2 {dx 2 + dy 2 ), 
with A = 4>{r)dt, 



9(f) 
%j) = ip( r )- 



(2) 



Its zero temperature limit has been studied in detail in 
[lH . Here, we will focus on the m 2 = case as the metric 
is non-sigular near the horizon (r = 0): 
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where ■0o,V ; i:Xi and ffi are all functions of a and xoi 
and their explicit forms are omitted here. Note that the 
metric becomes an AdS4 near the horizon. At the infinite 
boundary (r — * oo), 
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ip = J + 



(4) 



where p and p have the holographic interpretation as 
chemical potential and density, while J and (Ob) are the 
source and the condensate of the bosonic operator Ob 
which is dual to tp and living on the boundary. We use 
the shooting method to find the values of a and \o that 
gives J = (sourceless condition) and a chosen value of p 
(or p) at the boundary. In the above choice of boundary 
condition on ifi we have set the conformal dimension of 
Ob as As = 3, while g and x asymptote to an AdS4 at 
r — > oo . 

Keeping p constant, we obtain the following numerical 
results by fitting, 



a = 2.62, /^(l 

V i 



,-9/0.8 



)• 



(O b ) = 0.025/xV tanh(<? - 1.05) 



(5) 



The p dependence follows directly from dimensional anal- 
ysis. For constant p, p is monotonically decreasing in q 



while (Ob) is monotonically increasing in q. From this 
we know that the interaction is attractive. Because it is 
increasingly difficult to approach small q using the shoot- 
ing method, our fit cannot determine precisely the critical 
value of q below which the condensate vanishes nor the 
order of this phase transition. 

Fermionic Spectral Function- The fermionic spectral 
function is a powerful tool to probe the properties of the 
superconductor. We study the two point correlator (a 
retarded Green's function Gij) of the fermionic operator 
Of, and we will focus on the spin averaged spectral func- 
tion, which is defined as the imaginary part of the trace 
of Gr (in the space of the energy oj and the momentum 
k) 



A(uj, k) = Im(TrGn) 



(6) 



Note that Gr can either be computed using the CFT or 
through the fermionic field W which lives in the bulk but 
couples to Of at the boundary. In the second approach, 
Of can turn into \& '. Hence, the two point function of Of 
is related to the two point function of \& at the bound- 
ary. One can show that summing the corresponding mo- 
mentum modes in the r direction gives A(u>, k) = 2 as 
| a; | — > oo. For comparison, a free fermion spectral func- 
tion on the boundary of a 4D flat space is 



A(u, k) = 2 M /a/w 2 - k 2 6 (M - |fc|) 



(7) 



In our case, the information needed to compute Gr 
can be obtained by solving the Dirac equation. With \E' 
minimally coupled to the gauge field A^ in the gravita- 
tional background of Eq.(J5]), the Dirac equation can be 
recast into the flow equation [l6| 



d r £± = -2mFy/g~ii£,± T (k T u(r)) ±(k± u(r))£±, 




with u(r) 



-9tt 



-(u + q e <j>{r)), 



(8) 



where q e is the fermion charge. The mass tuf of ^ is 
related to Ap, the mass dimension of Of, as Af = 
niF + 3/2. Because of the strong interaction nature of 
the problem and the lack of supersymmetry in our back- 
ground, it is not clear how the charge and mass of the 
scalar field is related to those of the fermion field. Never- 
theless, we will focus on a possible scenario: q = 2q e and 
nip = m = 0. Thus, the scalar field has the same charge 
as a fermion pair, and the dimension of the bosonic op- 
erator Ob is just twice of the fermionic operator Op 
(A B = 2A F ). With m F = 0, we have 



TrG R = (£++£- 



(9) 



In the calculation, \ could approach some nonzero constant at 
r —* oo, implying a gravitational redshift. One can rescalc t —* 
e — x/2j anc j ^ e x/2^, t Q £ ne Einstein frame where the speed 
light is normalized to 1. 



Results and Discussions- Figure 1 shows the spectral 
function A(u>, k) of (q, p) — (1.7, 10) for sampled values of 
k. A(u), k) is even in k as (—k) = £-(k) in Eq.© for 
mp = 0. So we only need to focus on positive k. We see 
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FIG. 1: A typical fermionic spectral function with the 
peak-dip-hump structure for k = 0.7 (black dotted) 0.9 (red 
dashed) 1.5 (blue line) 2.5 (purple dotdashed). A sharp peak 
develops when k increases and the peak disappears above the 
Fermi momentum kp. Here, p — 10 and q — 1.7. 
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FIG. 2: 3D plot of the fermionic spectral function with the 
peak-dip-hump structure. Here, p = 10 and q — 1.7. 



at large \oj\ A(ui, k) asymptotes to 2, as mentioned above. 
There is an A(w, k) = region bounded by \ui\ < vq \k\, 
where vq is k independent and vq < 1. Outside this 
region, there appears a peak-dip-hump structure in the 
negative to (hole excitation) side when k is below some 
critical value kp. As k increases (but still below kp), 
the peak gets sharper and taller. Furthermore, both the 
peak and dip move toward smaller \lo\ while the hump 
moves toward bigger |a>|. The corresponding 3D plot of 
\ogA(ui, k) is shown in Fig. 2. It can be seen clearly 
that the peak becomes sharpest and tallest close to kp 
and disappears just above kp. It is encouraging that 
this behavior is very similar to what is observed in the 
superconducting phase of several HTSC materials. 

We identify kp as the Fermi momentum, because closer 
to the Fermi surface a hole excitation has fewer neigh- 
boring particles to annihilate with, so that its peak is 
sharper. We also identify the excitation energy of the 
sharpest peak as the gap A. It is consistent with the 



conventional definition of the lightest quasiparticle exci- 
tation energy. Numerically, for a wide range of q and 
fi (or p) values, we find that the gap (in the fermionic 
spectrum) is related to the bosonic condensate as 

0.6 



A 



P 



(10) 



This A oc (Ob) relation is also found in the BCS theory. 
It shows that the dynamics of the fermion and boson 
fields in the holographic theory are intimately related. 

It is instructive to show parametrically how the above 
relation could emerge out of the bulk theory. In the bulk, 
the term |qvl?/;| 2 in Eq.([T]) gives the gauge field a mass 
q 2 |i/'| 2 (Higgs mechanism). The exchange of one mas- 
sive gauge boson gives a localized potential between two 
fcrmions which can then be written as an effective four 
fermion contact interaction with coupling C ~ l/IV'l 2 - 
This interaction can contribute to the fermion self energy 
diagram with the fermion mass (i.e. the gap A) deter- 
mined self-consistently. This gap equation is solved when 
the fermion loop diagram / = 1/C. Since the fermion 
loop diagram has two fermion propagators and a four 
momentum integral, dimensional analysis gives / oc A 2 . 
Thus, we have A oc \ip\ oc (Ob)- 

Now we reach an interesting conclusion from the above 
discussions. From the analogy to the BCS theory, the 
bosonic condensate is the consequence of the fermion 
pairing. On the other hand, the condensate is also caus- 
ing the pairing, because the condensate gives the gauge 
boson mass through the Higgs mechanism which is im- 
portant for the fermion pairing. 

To understand the A(u), k) — region bounded by 
\uj\ < vq \k\, it is easier to go to the q — > oo limit. By 
rescaling A — > A/q and tp — > ip/q in Eq.|T]), A and tj) 
decouple from the Einstein equation; their back reaction 
can be neglected so the background is just AdS4 [Tl |. 
In this "probe limit," the bulk fermion satisfies the 4D 
Klein-Gordon equation 

[-V 2 + fc 2 -( W + ge 0(r)) 2 l 



ip = 0. 
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This is similar to a Shrodinger equation with a potential 
V(r) = —(us + q e cf)(r)) 2 and energy —k 2 . (j>(r) is zero at 
r = 0, and it increases smoothly and monotonically to 
fi at r — > oo. One can show that when \uj\ < \k\, tp is 
a bound state in the r-direction. This implies the bulk 
fermion propagation is not on-shell and hence Im(G/?) 
and A(lo, k) vanish. Away from the probe limit, the back 
reaction modifies the metric near r = but not at large 
r. Near r — 0, the light velocity factor v(r) = \/—gu/gtt 
starts to deviate from unity and decreases with decreas- 
ing q. The potential of the Schrodinger equation is then 
modified to V(r) — —v 2 (r)(u + q e (j)(r)) 2 , which makes 
the A{uj,k) — region to be bounded by \w\ < v \k\ 
with vo — v(0) (vq = 1 in the probe limit). In Figs. 1 
and 2, v = 0.15 for q = 1.7. 

Now we discuss the physics of the hump. The hump 
has the dispersion relation 



^hump 



(12) 
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and its width is ~ 0(q e p). The dispersion comes from 
the fermion gauge coupling to A t on the boundary, which 
gives an energy shift u> — > lo + q e p to the dispersion lj = 
—k of massless hole fermions. This is a mode of fermion 
propagation parallel to the boundary. More specifically, 
near the boundary, the bulk fermion satisfies Eq. (fTTj) . 
Thus, if <f>(r) were a constant p, the spectral function 
would be similar to Eq. ([7]) with a sharp peak at w = 
—q e p—k. In reality, 4>{r) almost remains constant except 
at small r. This curvature introduces a smearing effect 
of the sharp peak and sets the scale of the hump width. 
The physics picture is that the gradient of cf>(r) shows 
that the electric field localizes near r — 0. This is just in 
analogy with the Meissner effect which is a consequence 
of the Higgs mechanism. The photon mass q\ip\ sets the 
scale of the hump width. Numerically, this scale is ~ qp 
(p is the only dynamical scale in the problem) so the 
width is as large as LOhump at smaller k. This is indeed 
what we have seen in Fig. 1. 

In the above example, we have seen that while the size 
of p does not affect the physical properties we discussed 
above, the size of q could cause more significant changes 
by controlling the geometry. Another example is how 
the peak position changes with respect to k. In Fig. 1, 
the peak moves toward smaller \w\ with increasing k. In 
general, this is the case for q < 2.6. For q > 2.6, the peak 
moves toward bigger |u>| when k increases. This might 



imply that the system could go through a quantum phase 
transition by changing q, which requires further studies. 

In conclusion, we have studied the spectrum of 
fcrmionic excitations in a gravitational model which has 
taken account of the back reaction of the condensate. 
The system exhibits the PDH lineshape in the fermionic 
spectral function which shows great similarity to that in 
the ARPES experiments of HTSC. In addition, we have 
found a linear relation between the gap in the fermionic 
spectrum and the condensate, indicating the condensate 
is formed by fermion pairing. This is highly non-trivial, 
and it opens possibilities to study new microscopic mech- 
anisms for superconductivity without explicit pairing in- 
teractions. 
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